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INTRODUCTION 


An  exact  involute  surface  (EIS)  is  defined  by  the  relation  [1] 

r  sin  a  =  c  (1) 

where  r  is  the  radial  coordinate,  a  is  the  arc  angle,  and  c  is  called 
the  involute  constant.  Angular  coordinate  6  is  given  by 

A 

0  -  =  —  (2  -  z  )  +  a  +  cot  a  -  a  -  cot  a  (2) 

o  c  o  o  o 

where  the  point  defined  by  r  =  rQ  (a  =  aQ) ,  0  =  0Q,  z  =  zq  is 
simply  a  base  point  on  the  EIS  and  A  is  a  constant  related  to  the 
number  of  plies  N,  ply  thickness  t,  and  constant  c  by 


■  «S  •  •)' 


(3) 


The  positive  directions  of  the  various  geometric  variables  are 
shown  in  Figure  1*  where  it  should  be  noted  that  the  (0,r,z)  axes 
form  a  right-handed  triad. 

Other  important  geometric  parameters  of  the  EIS  are  the  tilt 
angle  y,  surface  angle  p,  and  helical  angle  <f>,  which  are  shown 
in  their  positive  directions  in  Figure  2.  These  are  defined  by 


dr 


tan  r  -  ii 


cosa 


(4) 


cos  ip  =  -  sin  a  sin  y 


(5) 


and 


<p  =  <p  +  <p  -  <p  +  — — r —  [  £  (z-z  )  +  cot  a  -  cota  ]  (6) 
°  (A^+l)  **  co  o 


where  <p  represents  the  angle  from  the  local  meridional  direction  to 
a  line  in  the  EIS  that  was  originally  straight  (in  the  ply  pattern) , 
such  as  the  warp  direction,  and  <t>Q  is  the  corresponding  angle  at 
the  base  point. 

^Figures  are  located  at  end  of  report. 
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The  developed  configuration  (  or  ply  pattern)  of  the  EIS 
is  defined  by  the  parametric  relations 

tan  X  =  |  (^±i2)  ^ 


r  -c 


A2z 


9  -  - , - r7_  -  X  -  cot  X 

c|X+l)  ' 

R  =  — 2  f (A2+l) (A2r2+c2) ]  h 

A * 


(7) 

(8) 

(9) 


where  the  ply  pattern  coordinates  R,  6  are  shown  in  Figure  3. 

The  foregoing  equations,  (1)  -  (9),  completely  define  the 
configuration  of  the  EIS  and  its  ply  pattern,  and  require  only 
the  prescription  of  constants  c,  N,  t  (or  equivalent  information) , 
initial  values  0  ,  ,  and  the  r,  z  profile  of  the  body  of  revolution 

being  generated  by  the  involute  method.  The  body  of  revolution 
itself  need  only  possess  piecewise-continuous  boundaries.  Equations 
(1)  -  (9)  guarantee  that  the  volume  of  the  body  will  be  completely 
and  perfectly  filled  by  layers  of  ideally  flexible  material  having 
constant  thickness  (aside  from  highly-localized  boundary  regions 
where  the  thickness  varies  to  avoid  "steps").  However,  exact 
involute  construction  is  not  widely  used  in  the  manufacture  of 
practical  components,  such  as  rocket  nozzle  exit  cones.  Rather, 
various  ply  pattern  models  based  upon  analytical  and  empirical 
methods  are  employed,  including  interleaving,  which  consists  of 
the  use  of  plies  of  (usually)  two  different  shapes  and  thicknesses 
in  the  same  body.  Thus,  it  is  important  to  mathematically  describe 
the  geometry  and  to  define  relations  that  measure  the  degree  of 
error  induced  via  these  approximate  methods,  which  is  the  purpose 
of  the  present  communication. 
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A  comment  on  the  organization  of  this  report  is  in  order  at 
this  time.  In  the  next  section,  we  shall  present  the  equations  that 
define  the  involute  surface,  i.e.,  the  three-dimensional  surface 
formed  by  a  ply  within  the  body  of  revolution,  and  the  associated 
ply  pattern  according  to  contemporary  procedure  [2,  3,  4],  This 
procedure  consists  of  the  analysis  of  involute  surface  "strips"  of 
finite  dimension  and  involves  iterative  solutions  of  the  basic 
equations.  Some  new  results  for  the  description  of  the  helical 
angle  distribution  and  ply  pattern  analysis  in  the  presence  of 
interleaving  are  also  given  there.  New  work  that  leads  to  precise 
(closed  form)  solutions  of  the  basic  equations  is  presented  in  a 
subsequent  section.  This  is  done  in  the  interest  of  clarity,  and 
also  to  establish  a  complete  formulation  by  the  finite  strip 
method  since  the  ply  pattern  analysis  does  not  lead  to  a 
convenient  closed  representation. 
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THE  START-LINE  METHOD 


a)  INVOLUTE  SURFACE  STRIP 

The  model  employed  herein  assumes  that  the  involute  surface 
consists  of  strips,  each  of  which  obeys  eqs.  (1)  -  (9)  with 
piecewise  constant  c,  such  that  consecutive  strips  are  continuous 
at  a  single  point.  The  consecutive  strips  possess  a  common 
terminal  plane  z  =  constant.  This  model  leads  to  perfect  filling 
of  volume  if  the  strip  dimension  is  infinitesimal  but  results  in 
discontinuities  in  the  form  of  "plateaus"  in  the  involute  surface. 
Furthermore,  we  introduce  the  concept  of  a  start- line,  as  described 
by  Savage  [3]  and  Stanton  [4]  to  control  the  local  orientation 

of  the  involute  surface.  The  start-line  is  simply  a  meridian 

(0  =  constant)  of  the  involute  surface  and  is  only  required  to 
be  piecewise-continuous.  It  seems  that  most,  if  not  all,  of 
the  ply  pattern  design  approaches  used  in  practice  can  be  derived 
by  appropriate  tailoring  of  the  start-line  [4].  Furthermore,  we 
shall  use  the  term  hypothetical  involute  surface  to  refer  to  the 
involute  surface  conceived  in  the  start-line  method,  which  cannot 
be  formed  from  sheets  of  continuous  material  because  of  the 
discontinuities  required.  It  will  be  necessary  to  refer  to  four 

specific  values  of  radius;  namely,  inner  and  outer  radii  r^.  and 

rQ,  respectively,  start-line  radius  rg,  and  for  the  case  where 
interleaving  is  present,  rc  is  the  intersection  of  two  regions 
having  different  values  of  N  and  t.  Each  of  the  four  radii  are 
functions  that  depend  on  axial  coordinate  z.  Since  the  constants 
N  and  t  usually  appear  in  the  form  of  their  product,  we  let 

M  =  Nt  (10) 
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We  shall  consider  a  strip  of  the  involute  surface  as  the  region 
formed  by  z^  <  z  <  z^.  Parameters  will  carry  a  subscript  q 
(q  =  1,2)  that  corresponds  to  thickness  parameter  M  if  inter¬ 
leaving  is  present.  Our  convention^-  is  that  q  =  1  corresponds  to 
the  region  r  >  rc  if  rg  >  rc  and  r  ^  rc  otherwise.  We  cannot 
permit  the  quantity  rg  -  rc  to  change  sign  in  a  given  problem, 
although  it  may  vanish.  Also,  a  subscript  G(G  =  1,0, S,C)  will  be 
used  to  refer  to  any  of  the  four  special  values  of  r  discussed 
above.  If  the  value  of  r  is  arbitrary,  no  G  subscript  will  be 
displayed.  It  will  also  be  understood  that,  unless  otherwise 
noted,  the  following  equations  are  only  valid  within  a  single 
strip,  hence  no  index  for  strip  identification  will  be  employed 
unless  needed  for  clarity.  Finally,  we  introduce  the  representation 

PGq  (VZk>  *  aGq  <zj>  +  TOt  »Gq  <V  ’  “Gq  (zk>  ‘  COt  °Gq  ,zk>  1111 
and  we  assume  that  rg(Z2)  -  rg(z^)  i®  not  negative. 

We  now  let  the  strip  geometry  be  defined  by  the  governing 
involute  equation 

r  sin  a  =  c  (12) 

q  q 


where  c^  is  piecewise-constant.  While  r,  0,  and  z  are  all 
continuous  at  r  =  rc,  arc  angle  a  is  discontinuous  there, 
let 


4*2c  2 


Aq=  <' 


3-  -  1) 


We  also 

(13) 


The  constants  cannot  be  arbitrarily  prescribed  when  using  the 
start-line  approach.  Rather,  c ^  is  computed  such  that  it  satisfies 


FS1  (zl'  z< 


(14) 


^Interleaf  terminations  on  planes  z  =  const,  are  not  treated  in 
this  work. 
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in  conjunction  with  (12)  and  (13)  with  q  =  1.  These  relations 
guarantee  that  the  points  on  the  hypothetical  involute  surface 
corresponding  to  r  =  rg  at  z  =  z^ ,  z^  lie  on  the  same  meridian  9 
=  0  .  A  good  first  approximation  for  the  solution  of  this  system 
of  equations  is  given  by 


where 


2 

C1  a 


PG  = 


1  +  Pc 


4*2/Mi2+ps2/rs2 


rG(22)~rG(zl) 
Z  2  "  Z1 


(15) 


(16) 


and  the  subscript  G  has  the  same  meaning  described  earlier. 

The  constants  are  governed  by  continuity  of  the  hypothetical 
involute  surface  at  r  =  rc,  z  =  z^,  z^r  which  lead  to 

A  A 

FC2(zl'Z2)  ”  FC1(Z1'Z2*  =  (cj  "  ^  (Z1  "  Z2J  (17) 

along  with  (12)  and  (13).  We  may  use 

M2 

C2  *  C1  (18) 

as  a  good  first  approximation  for  the  solution  of  (17),  (12),  and 

(13) ,  although  solutions  do  not  exist  for  arbitrary  values  of  the 

input  parameters.  The  involute  surface  strip  is  thus  defined  by 
A 

0  (z)  =  9g  +  -S'  (z-Zj^)  +  asl(zl)+  cot  oig^z^-oi  (z)  -  cot  a  (z) 


+  [°C2(*i)  +  cot  «C2(zi)  "  “ci^l^  “  cot  «C1<Z1)]  S  q2  (19) 

which  is  continuous  at  z  *  z z^.  Here,  6 is  the  Kronecker 
delta;  specifically,  5  22  =  1  and  5 12  “  0 ’  The  tilfc  and  surface 
angles,  defined  .  _ 
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and 


tan  y 

q 


A 

_a 

cosa 

q 


(20) 


cos  <j> q~~  sin  “q  sin  (21) 

like  a  ,  contain  discontinuities  at  r  =  r.,. 

We  next  consider  the  distribution  of  helical  angle  <j>,  where 
we  restrict  our  attention  to  the  case  in  which  the  warp  (and  fill) 
fibers  of  the  regions  q  =  1,  2  in  the  ply  pattern  are  parallel. 
Within  each  region  of  the  hypothetical  involute  surface,  the 
helical  angle  is  given  by  an  equation  of  the  form (6) ,  however,  a 
discontinuity  exists  at  r  =  rc.  To  define  the  magnitude  of  the 
discontinuity,  we  consider  a  strip  of  infinitesimal  height  such 
that  z^  =  z  and  z^=  z  +  dz.  In  this  strip,  a  crease1  occurs  at 
r  =  rc>  As  one  may  observe  by  creasing  a  piece  of  lined  paper, 
the  angles  between  the  warp  direction  and  crease  line  on 

.  .  ,  ❖  A  iA 

either  side  of  the  crease  line  are  equivalent.  Letting  (i,  j,  k) 
represent  unit  vectors  along  the  (0,  r,  z)  directions,  respectively, 
at  rc,  the  infinitesimal  vector  along  the  crease  line  is  given  by 

d£  =  -rd0i  +  dr]  +  dzk  (22) 

But,  differentiation  of  (19)  leads  to 


which. 


on  use  of 


A  2 

d0  =  -3.  dz  +  cot* 
c 

q 

(12),  becomes 


Cq 


(23) 


c  d0  =  A  dz  -  cos  ar  dr  (24) 

q  q  <~q 

Now  the  unit  vectors  along  the  warp  direction  on  either  side  of 
rc  are  given  by  wq  where  (see  eqs.  8,  12,  13  of  Reference  [2]) 

1  In  actual  involute  bodie  ,  the  crease  appears  to  be  smoothed  out. 
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ttur  -A- 


w  —  - 

q  sinij' 


Cq 


tcos  aCq  Sin  ^Cq  1  +  (8in  YCq  COS  *Cq  sin  ^Cq 


-  sin  aCq  sin  <»Cq  cos  YCq)  j  +  cos  YCg(cos  4>Cqsin  *Cq 


(25) 


-cos  ^Cq  sin  <|>Cq)£3 


Letting  =  ft2,  or 


d£  •  =  d£  •  w2 


(26) 


where  the  dot  signifies  a  scalar  product,  leads  to  the  relation 

f2  sin  <)>c2  +  g2  cos  4>c2  =  f1  sin  <j»cl  +  gx  cos  4>cl  (27) 

where 


1  2 

f  =  —7 — : [  (cosa„  cota„  -  sina„  cos  y„  )  p  -A  cota_ 

q  sinijj  Cq  Cq  Cq  'Cq  q  Cq 

eg 


-  cos^Cq  cos^Cq1 


9q  *  cos^Cq  +  PC  sinYCq 


(28) 


(29) 


recalling  that  pc  is  given  by  (16).  Thus,  4>c2  is  defined  by  the 
solution  of  (27)  ,  once  4>c^  has  been  established.  The  solution 
of  (27)  is  given  by 

_!  Q, 

(30) 


-1  ®l 

♦C2  =  U2i  COS  [Qj  COS  {*Cl  ‘  U,l)] 


where  the  positive  sign  is  chosen  when  >_  co^,  otherwise  the 

negative  sign  is  used.  In  (30),  we  have  made  the  replacements 

f 

tan  u  =  -2  (31) 

q  gq 


Q  =  (f  2  +  g  2)is 

q  q  yq  ' 

and  we  have  assumed  that  0  <  <}>-  -  ui  <  n. 

-  Cq  q  - 


(32) 
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In  the  derivation  of  (30),  we  have  assumed  the  strip  has 
infinitesimal  height  Z2  -  z^  =  dz.  For  the  case  of  finite  height, 
we  use  the  same  relations  as  an  approximation.  Furthermore,  we 
choose  to  write  (30)  at  each  point  r^,  ( z)  in  order  to  develop 
the  general  relation  for  4>2(z)  ,  although  other  interpretations 
for  this  approximation  are  possible,  such  as  writing  (30)  only 


at  rc(z^)  and  then  utilizing  an  equation  of  the  form  (6)  for  q  =  2 
with  base  point  (rQ,  zq)  replaced  by  (rc(z^),  z^).  The  general 


expressions  for  <J>  thus  become 


4»l(z)  =  *sl(z1)  +  ♦siUi)  -  ^(z)  +  -—2— t—  (z-zx)  + 

(A^  +1)  2  1 


cotagl(z1)  -  cot  a1(z)] 


1>2(z)  =  <*>C2(z)  +  ^C2(z^  ~  *2*Z*  + 


cot  (z ) ]  which  follow  from  (6) 


- - - -  [cot  a  (z)  - 

(A22+l)!5  C2 


b)  PLY  PATTERN  OF  A  STRIP 

We  now  turn  our  attention  to  the  ply  pattern,  or  developed 
view,  of  a  basic  strip  of  th6  hypothetical  involute  surface.  For 
this  purpose,  we  employ  a  set  of  equations  of  the  form  (7)- (9)  to 
represent  each  of  the  mapped  regions  q  =  1,2.  In  other  words, 
consider  two  systems  of  plane  cylindrical  coordinates  (R^,  0  ) , 
along  with  associated  parameters  A  .  If  the  involute  construction 
procedure  were  exact,  we  would  position  the  origins  and  orient  the 
axes  of  the  two  coordinate  systems  such  that  the  ply  pattern  would 
be  continuous  along  the  crease  line.  But  all  solutions  other  than 
an  exact  involute  surface  are  approximate,  hence  this  continuity 
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condition  cannot  be  perfectly  satisfied.  Our  approach  to 
approximate  continuity  consists  of  selecting  the  two  R,  0 
coordinate  systems  such  that  the  ply  pattern  is  continuous  at  the 
mapped  point  corresponding  to  rc(z^)  while  the  distance  between 
the  two  mapped  images  of  r^(z2)  is  minimized.  In  this  manner, 
three  parameters,  corresponding  to  relative  translation  and 
rotation  of  the  two  coordinate  systems,  can  be  evaluated.  We 
shall  subsequently  employ  a  similar  algorithm  to  treat  strip-to- 
strip  continuity  of  the  ply  pattern.  Thus,  we  shall  use  subscripts 
a,  b  in  lieu  of  1,  2  in  order  to  avoid  specific  dependence  on  the 
region  index  q  in  what  follows. 

The  two  plane  cylindrical  coordinate  systems,  (R  ,  0  )  and 
(Rb#  /  are  shown  in  Figure  4(a),  as  well  as  a  cartesian 
system  (x,y) .  For  our  purposes,  we  shall  regard  (x,y)  as  being 
fixed  in  space,  while  the  former  two  vary  in  orientation  and 
position  according  to  the  strip  and  region  under  consideration. 

From  Figure  4(a),  we  observe  that 

x  =  xa  +  Ra  COS  (  Ga  '  V  =  xb  +  *b  cos  (  %  ~  V 

(34) 


y  =  ya  +  Ra  sin  ( 0a  -  ©a>  =  Yb  +  %  sin  (^  -  y 


for  a  point  P  in  the  ply  pattern  which  has  coincident  images 
in  the  (Ra,  ©a)  and  (R^,  0^)  coordinate  systems.  If  the  corresponding 
images  of  a  mapped  point  do  not  coincide,  they  are  separated  by 
a  distance  L  which  is  given  by 


h2  =  lxb-xa  +  Rj/cos  (0*b-%)  ‘  Ra*  cos  (0a*  "  ®a)]2  +  ^b'^a  + 
"b*  sin  (^*  -  \)  -  Ra*  sin  ( ©a*  -  0a)]2 
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It  should  be  noted  that  eq.  (37)  always  has  two  roots.  The  correct 
root  is  the  one  which  minimizes  L  in  eq.  (35). 

Thus,  to  satisfy  the  aforementioned  continuity  conditions  for 
adjoining  regions  q  =  1,2,  we  simply  let  a  =  1,  b  =  2  (or  vice 
versa)  and  make  the  replacements 


Rq  =  RCq(zl)'  ®q  "  GOq(zl) 


(z  )  , 
2 


in  (38),  which  in  turn  is  substituted  into  (37)  and  thence  into 
(34),  so  that  the  coordinates  of  the  ply  pattern  strip  become 
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(40) 


x  =  x  +  R  cos  ( 0  -  ©  ) 

3  <3 

y  =  y  +  R  Sin  ( 0  -  ©  ) 

M  q 

according  as  the  image  point  (R,  0)  corresponds  to  q  =  1  or  q  =  2. 

In  (40)  ,  we  let  (R,  0  )  represent  the  ply  pattern  parameters  of 

an  arbitrary  point,  which  may  lie  in  either  region. 

c)  SPECIAL  CASE:  CYLINDRICAL  INVOLUTE  STRIP 

In  the  usual  case,  R  and  0  are  given  by  eqs.  (7)  -  (9).  An 
exceptional  case  occurs  however,  when  the  hypothetical  involute 
surface  strip  is  cylindrical,  i.e.,  y  =  A  =  0.  In  this  case, 
eq.  (9)  is  not  valid  and  the  ply  pattern  degenerates  into  a 
trapezoid  in  which  the  mapped  edges  corresponding  to  the  planes 
z  =  constant  are  parallel.  Furthermore,  the  length  of  each  of 
these  mapped  edges  is  equal  to  its  arc  length  in  the  involute 
surface.  Thus,  referring  to  Figure  4(b),  we  redefine  the  quantities 
R  and  0  for  this  special  case  by 

2  ,  [r2 (z) -rr2 (z. ) ] 2 

RCYL(z)  =  (z-zl}  +  - —2—1 -  (41) 

4c 

q 

sin  W*1  -  £a>  (42) 

where  G  may  be  set  equal  to  S  if  q  =  1  and  C  if  q  =  2  for 
convenience.  In  this  way,  the  origin  of  R,  ©  coordinates  in  Figure 
4(b)  is  placed  at  the  image  of  rg(z^)  if  q  =  1  and  r^z^)  if  q  =  2. 
With  the  interpretation  of  (41)  and  (42)  for  any  cylindrical 
involute  region,  the  previous  (and  subsequent)  ply  pattern  equations 
are  all  valid. 
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d)  CONNECTIVITY  OF  STRIPS 

The  preceeding  equations  (10)  -  (21)  define  the  hypothetical 
involute  surface,  or  spatial  configuration  of  a  ply  within  the 
generated  body  of  revolution  conceived  in  the  start-line  method. 

As  mentioned  earlier,  the  strips  are  only  connected  along  the 
start-line  itself.  Internal  continuity  conditions  between  regions 
q  =  1,2  have  already  been  established.  Thus,  for  piecewise- 
continuous  functions  rG(z),  these  equations  represent  a  complete 
system.  It  is  only  necessary  to  recall  that,  for  consecutive 
strips,  the  previous  value  of  z becomes  the  new  value  of  z^ffor 
the  next  strip) . 

In  order  to  define  the  helical  angle  distribution  via  eqs. 
(33),  an  additional  constraint  is  required,  for  example 


(i) 


(i+1) 


<(’S1(Z2)  =  *Sl(zl> 


(43) 


where  we  have  introduced  a  strip  index  as  a  superscript  since 
the  relation  involves  consecutive  strips.  However,  because  of 
the  inherent  approximation  of  the  start-line  method,  (43)  must  be 
viewed  as  an  assumption.  It  may  be  just  as  accurate  to  assume 
that  4>sl(z)  is  a  constant.  For  typical  exit  cones  in  which  a  is 
quite  small,  the  difference  between  these  two  assumptions  would 
be  negligible. 

Unlike  the  hypothetical  involute  surface,  however,  the  ply 
pattern  equations  (34)  -  (42)  do  require  explicit  connectivity 
relations  to  define  strip-to-strip  continuity.  Here,  the  above 
algorithm  involving  eqs.  (38),  (37),  and  (34)  is  repeated,  however, 
in  place  of  (39)  we  use 
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(i)  (i) 

Ra  =  RS(z2)'  °a  =  es(z2} 


(i+1)  (i+1) 

Rb  =  RS  ^ z  1 )  '  ^3  ®s^Zl^ 


*  (i)  *  (i) 

Ra  =  RG(z2}'  0a  "  °G(Z2)  (44) 


*  (i+1)  *  (i+1) 

%  =  W'  %  =  W 


where  G  *  S  may  be  chosen  arbitrarily.  Convenient  choices  are 
G=CorG=  IorO,  corresponding  to  the  farthest  point  from 
the  start-line.  The  ply  pattern  coordinates  are  thus  given  by 


(i)  (i)  (i)  (i)  (i) 

x  =  x  =  x  +  R  cos  (  ©  -  0  ) 


(i)  (i)  (i)  (i)  (i) 

y  s  y  =  yq  +  R  sin  {  0  -  0g) 

(i)  (i) 

in  a  common  coordinate  system  (x,y).  The  values  of  x  ,  y  ,  and 

(i)  q  q 

(j  for  single  fixed  values  of  i  and  q  may  be  chosen  arbitrarily. 

The  framework  provided  by  eqs.  (10)  -  (45)  now  constitutes  a 
complete  system  of  recursion  relations  such  that  of  the  sub¬ 
sequent  strip  corresponds  to  z of  the  present  strip.  Input  data 
for  this  model  requires  prescription  of  the  z  -  dependent  radii, 
rJf  rQ,  rc,  and  rg  (usually  given  in  digital  fashion),  M^lq  =  1,2), 
the  start-line  meridian  9g,  and  the  helical  angle  0  at  a  single 
point  on  the  start- line.  In  the  event  that  no  interleaving  is 
present,  one  may  simply  neglect  equations  related  to  the  region  q 
=  2  or  set  M2  =  in  the  general  equations.  The  start-line 
procedure  always  yields  discontinuities  in  the  ply  pattern  and  the 
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involute  surface  unless  an  exact  involute  surface  is  employed. 

In  the  latter  case,  the  start-line  is  no  longer  arbitrary  and  it 
is  only  necessary  to  use  eqs.  (1)  -  (9).  If  the  exact  involute 
surface  is  also  cylindrical,  then  eqs.  (41)  and  (42)  can  be  used 
in  place  of  (7)  -  (9) . 


EXPLICIT  RELATIONS  FOR  THE  START-LINE  METHOD 


In  this  section,  we  shall  demonstrate  that  the  equations 
defining  involute  parameters  c^,  i.e.,  eqs.  (12)  -  (14)  and  (17), 
can  be  solved  in  explicit  (rather  than  iterative)  fashion  provided 
that  we  model  the  involute  surface  by  strips  of  infinitesimal 
dimension  -  z^. 

Letting  z^  =  z  and  =  z  +  dz ,  we  observe  that  (14)  may  be 
expressed  as 

A1 

Fsl(z,  z  +  dz)  +  ^  dz  =  0  (46) 


while  (11)  gives 


Fgi(z,  z  +  dz)  =  cot  asldotsl 


(47) 


where  agl  is  the  arc  angle  at  rg(z).  Now  substituting  (47)  into 


(46)  after  use  of  (12)  and  its  derivative  leads  to 

dr„ 


2  _  2V  I* 


A^dz  = 


<rS  C1  > 


(48) 


whence,  substituting  (13)  into  (48)  and  solving  for  c^,  we  get 

o  1  +  m„2 


'1  .  2  .  2  2.2 
4tt  /m-j^  +  mg  /rg 


(49) 


where 


dr 


mG  = 


_G 

dz 


(50) 


We  also  see  that  mQ  is  equivalent  to  tan  yg  by  (4),  however,  the 
form  (50)  can  be  expressed  in  terms  of  given  information,  thus  is 
more  convenient  for  use  in  (49).  Equation  (49),  which  has  been 
given  earlier  as  a  good  first  approximation  for  finite  strip  height 
shows  that  c^  in  a  practical  exit  cone  built  by  the  start- 
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line  method  never  differs  greatly  from  its  value  for  cylindrical 
M1 

involutes  . 

In  similar  fashion,  we  may  derive  an  explicit  expression  for 

c_.  In  this  case,  (17)  becomes 
^  A  A 

F  (z ,  z+dz)  -  Fri(z,z+dz)  =*  (-^  -  — )  dz  (5) 

C2 

Using  (11)  and  (12)  along  with  its  derivative  again,  we  find  that 


(rc2-°i2)i4 


.  2  2 . 
(rC  ~C2  ] 


2  2 

m  A  4  it  c_ 

]  ^  ( - y~  -  1) 

rC  C1  C2 


Solving  for  c2 ,  we  get 


where 


a.+a0  ^ 

c2  *  '-it1’ 


ax  =  (l+mc2)  (4ir2/M22+mc2/rc2)  +  (mc2-l)a42 

a2  =  2mca4  [  (l+rac2)  (4Tr2/M22-l/rc2)  -  el*)*  (54 

a3  =  (4*2/M22-mc2/rc2-a42)2  +  16Ti2mc2/M22rc2 

a,  =  ±  [(41t2c12/M12  -  1)^  -  (r^V  mc/rc) 

The  sign  ambiguity  in  (53)  is  caused  by  an  extraneous  root 
resulting  from  the  solution  procedure.  The  correct  root  is  the 
one  that  satisfies  eq.  (52) .  In  order  for  a  solution  of  (52)  to 
exist,  the  input  parameters  must  satisfy  the  relation 


1  a4 

M_  ^  2  TT  ( - y  +  - y) 

2  rc  i+”c 
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which  simply  requires  to  be  real.  Also,  if  the  region  q  =  1 
is  a  cylindrical  involute  while  q  =  2  is  not,  the  data  must  be 
such  that  c^  > 

Although  the  equations  governing  the  ply  pattern  may  be 
solved  explicitly  as  above,  the  results  are  very  cumbersome  and 
inconvenient,  thus  they  will  not  be  derived  here. 

I 

I 


DISCONTINUITY  STRAINS 


An  involute  surface  is  conceived  to  be  built  of  strips  in 
the  start-line  approach.  This  hypothetical  involute  surface,  in 
general,  contains  discontinuities  in  the  form  of  "plateaus"  at 
interfaces  z  =  constant.  Since  the  procedure  is  approximate  it 
is  important  to  quantify  the  error  in  order  to  compare  ply  patterns 
and  to  define  regions  in  which  ply  pattern  deformations  tend  to 
become  excessive.  Our  proposed  method  of  evaluation  of  ply  pattern 
accuracy  depends  on  the  introduction  of  "discontinuity  strains," 
which  are  defined  as  the  components  of  an  applied  strain  field  that 
would  transform  the  hypothetical  involute  surface  into  a  continuous 
sheet. 

We  again  consider  a  hypothetical  involute  surface  generated 
by  infinitesimal  strips  and  its  traces  in  the  planes  z,  z  +  dz  = 
constant.  However,  at  z  +  dz,  two  traces  are  considered  on  the 
plateau.  These  correspond  to  involute  constants  c  and  c  +  dc 
and  are  distinct  except  on  the  start-line,  where  the  two  traces 
intersect.  In  this  way,  we  establish  an  algorithm  suggested  by 
Figure  5,  where  the  six  points  1-6  are  introduced.  Corresponding 
numerical  values  are  used  as  subscripts  to  define  coordinates. 
Points  1,  2,  3,  and  4  all  lie  on  a  continuous  surface  with 
involute  constant  c,  while  5  and  6  are  on  a  surface  defined  by 
involute  constant  c  +  dc.  The  discontinuity  strains  are  the 
(hypothetical)  strains  introduced  in  the  transformation  removing 
the  discontinuity,  i.e.,  moving  3  to  5  and  4  to  6. 

In  general,  the  discontinuity  strains  will  not  represent  the 
actual  strain  field  in  the  fabricated  body  since  their  presence 
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would  depend  upon  specific  (artificial)  constraints.  However,  we 
may  state  that  a  continuous  involute  surface  would  result  from  their 
application.  The  discontinuity  strains  are  a  direct  measure  of 
the  errors  involved  in  the  use  of  the  start-line  method,  and  they 
may  be  useful  in  delineating  potential  areas  of  wrinkling  and 
severe  fabric  distorsion.  Thus,  the  discontinuity  strains 
should  be  considered  as  a  guide  to  good  ply  pattern  design,  rather 
than  as  a  precise  measure  of  material  deformation. 

Recalling  eq.  (19),  we  can  observe  that 

2 

d0  =  cot  add  lc,\ 


c  2  ( S', 

dz  =  -  £  cot  ada 
A 

where  the  region  index  q  is  dropped  since  it  is  immaterial  for 
the  present  analysis.  Thus  from  (56),  (57),  and  (12),  we  arrive 


r  =  r  -  £l9_ 

2  1  cosa. 


r,  =  r.  + 

3  1  cosct. 


_  r  cd9  +  Adz 

r4  rl  cosa^  cosa^ 


where  0.  =  6 =  0  and  0_=0.=0+d0. 

13  2  4 

It  is  now  necessary  to  define  the  location  of  points  5  and  6 
such  that  we  maintain  a  one-to-one  correspondence  between  them¬ 
selves  and  their  associated  points  3  and  4.  To  accomplish  this, 
we  observe  that,  despite  the  appearance  of  two  traces  at  z+dz, 
the  radii  at  the  end  points  of  the  two  traces  (rQ  and  rj)  are 
identical  since  rG  are  functions  of  z  only.  Thus,  the  ratio  of 
arc  lengths  along  the  two  traces  is  inversely  proportional  to  the 
ratio  of  their  involute  constants.  This  relationship  can  be 
preserved  by  taking 
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(59) 


r5  r  3 


r6  r4 


Thus,  the  transition  between  the  traces  at  z+dz  involves  a  change 
in  a,  but  not  r.  Using  (12)  then,  we  get 


da  _  tana 
dc  c 

at  z+dz,  while  at  the  same  level,  we  have 


d(a+cota)  dc  „  ,  dc  ,  2  2.Jj 

- dS -  =  "  “c"  cota  =  -  ~2  (r  -c  ) 

c 

Substituting  these  relations  into  (19)  leads  to 


where 


e5  =  e  +  <ae 


d6  =  [c  cota1  -  (rF2-c2)  **  ] 


and  rF  is  the  radius  at  the  point  where  the  traces  at  z+dz  intersect, 

which  is  r  if  q  =  1  and  r  if  q  =  2 .  We  also  get 
o  L 


eg  =  95  +  d0  (64) 

from  (60),  (61),  and  (19). 

To  define  the  discontinuity  strains,  we  begin  by  computing 
the  vectors  s^,  and  s,.  that  are  directed  from  point  1  to  the 
respective  points  as  shown  in  Figure  5.  In  terms  of  a  unit  triad 

AAA 

(i,  j,  k)  directed  from  the  origin  toward  point  1,  we  have 

s-  =  (r.  i  -  — f—  j )  d6 
2  1  cosa^ 


s  =  ( _ £L_ 

s3  'cosa, 


j  +  k)dz 


mm  S'  A  /v 

s-  =  r.dsi  +  j  +  dzk 

5  1  cosa^ 
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where  we  have  dropped  higher  order  terms  and  d3  is  given  by  (63). 
Thus,  the  strain  in  the  meridional  direction  (see  Figure  6)  is 
given  by 


e  =  —  -  1 
m  s3 


(66) 


which,  on  use  of  (65)  and  the  assumption  that  em  <<  1,  yields 


em  = 


(c  6'cota^) 

2 (A2+cos2a^) 


(67) 


where  the  prime  denotes  differentiation  with  respect  to  z.  We 
note  that  8  '  can  be  expressed  as 

8'  =  ^  [c  cotci.  -  (rp2-c2)55  ]  (68) 

c 

with  c'  given  by  differentiation  of  (49)  or  (53) ,  according  as 
q  =  1  or  2 ,  respectively. 

We  now  consider  the  strain  in  the  s-direction,  where  s 
corresponds  to  the  tangent  vector  in  a  plane  normal  to  the  z-axis. 
The  appropriate  deformation  involves  a  transformation  that 
carries  the  arc  connecting  points  3  and  4  into  that  connecting  5 
and  6.  The  ratio  of  these  two  arcs,  however,  is  the  inverse  of 
the  ratio  of  their  respective  involute  constants,  which  are  c  and 
c+dc.  Hence,  only  infinitesimal  strain  accompanies  this  trans¬ 
formation  and 


Eg  =  0  (69) 

For  our  third  deformation  measure,  we  shall  compute  the 
distortion  of  surface  angle  tp.  Letting  ft  represent  the  angle 
between  S2  and  s^,  we  get 


cos  ft 


s2  *  s5 
S2S5 


(70) 
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(71) 


which,  on  substitution  of  (65),  becomes 


cos  a  = 


(c  S'  cot  a^-AJsina^ 


[  (c  S'  cota^)  +A^+cos^a^]  ** 


Using  elementary  trigonemetric  relations  in  conjunction  with  (5) 
brings  us  to 


sinfil^-ft)  = 


| (A^+l) ^ (cS  1  cot^a^-A) +A [ (A+cS1  ) ^+1 ]  2| sina^cosa^  (72) 
(A^+cos^a^)  ^  [£  S  1  cota^)  ^4-A^+cos^a^]  ^ 


The  quantities  given  by  eqs.  (67),  (69),  and  (72)  may  be 
viewed  as  deformation  measures  in  a  non-orthogonal  coordinate 
system  since  the  angle  between  s  and  n  is  i|i,  rather  than  tt/2. 
These  quantities  may  be  transformed  into  discontinuity  strains  in 
the  orthogonal  system  (s,t)  shown  in  Figure  6  via  the  relati<'rs 

2 

Gm  =  etsin  'J'l+Y^sin^cosi^  (73) 


The  terms  involving  cot^  tend  to  be  small  in  practical  exit  cones 
since  \j>  is  only  slightly  larger  than  tt/2.  The  strains  (74)  may 
now  be  t ran i formed  into  arbitrary  coordinate  systems,  such  as  the 
warp/fill  axes,  by  use  of  the  standard  strain  transformation 
equations  of  linear  elasticity. 
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DEBULKING  KINEMATICS 


In  this  section  we  shall  derive  explicitly  the  equations 
referred  to  in  [1]  pertaining  to  the  analysis  of  debulking 
kinematics.  Specifically,  we  shall  consider  the  configuration 
of  the  involute  surface  corresponding  to  material  in  its  initial, 
undebulked  form  assuming  that  the  ply  pattern  based  upon  cured 
thickness  is  given,  starting  with  the  premise  that  it  is  desirable 
to  pressurize  against  tooling  that  controls  the  profile  of  either 
final  surface,  i.e.,  r^z)  or  rQ(z),  during  debulk.  The  treatment 
is  valid  specifically  for  exact  involute  construction.  Therefore, 
if  involute  construction  based  upon  the  start-line  method  is  employed, 
the  model  given  here  applies  strictly  to  a  strip  of  the  hypo¬ 
thetical  involute  surface. 

As  discussed  in  [1],  it  is  not  possible  to  design  a  ply 
pattern  for  an  arbitrary  body  of  revolution  such  that  one  of  its 
edges  will  remain  in  contact  with  a  control  surface  throughout 
a  debulking  procedure.  However,  approach  was  shown  that  permits 
an  edge  of  the  ply  pattern  to  lie  quite  close  to  the  control 
surface  during  debulking  over  a  significant  region  provided  that 
the  ply  pattern  contains  no  discontinuities  in  the  region.  In 
this  approach,  the  two  end  points  of  the  ply  pattern  boundary  in 
the  (smooth)  region,  which  are  termed  reference  points,  lie  on  the 
control  surface.  The  formulation  proceeds  as  follows: 

Determine  the  ply  pattern  corresponding  to  the  final  (cured) 
state  in  the  usual  manner.  We  let  the  initial  (undebulked)  position 


24 


of  point  1  coincide  with  its  final  location,  while  point  2  is 
merely  constrained  to  lie  on  the  control  surface.  The  coordinates 
and  other  parameters  evaluated  at  points  1  and  2  shall  be  denoted  by 
the  respective  subscripts . ^  The  usual  notation  is  employed  for 
the  cured  configuration,  while  the  same  symbols  with  asteriks 

will  be  used  to  refer  to  undebulked  configuration.  For  the  first 

*  *  * 

approximation,  let  the  initial  coordinates  (r2  ,  02  ,  z2  ) 
coincide  with  their  final  values.  Using  (19)  then,  we  arrive  at 

6  *  -  0.  =  -*  (z  *  -z.)  +  a  * +  cota  *  -  a  *  -  cota  *  (75) 

21c*  211  12  2 

for  a  first  approximation  of  c*.  Equation  (75)  can  be  solved  by 
iteration  starting  with  an  assumed  value 


Now  let  the  control  surface  be  defined  by 

r2*  -  r1  =  f  (z2*-z1)  (77) 

where  f  is  a  known  function.  Also  let  £  represent  the  distance 
between  the  images  of  points  1  and  2  in  the  ply  pattern,  or 

£  2  =  R12+R22-2R1R2cos  (©2-%)  (78) 

We  then  solve  for  r2*  and  z2*  by  iteration  using 

2 

(z  *-z  ) (A*) 

0*  =  — - - ^ - j—  -  X  *  -  cot  X  *  +  X  *  +  cot  X  *  (79) 

c*[(A*)  +1}H  2  21  i 

and 

(Rx*)2  +  (R2*)2  -  2R1*R2*cos©*  =  £2  (80) 

1  Similar  to  the  previous  section,  numerical  subscripts  are  only 
used  in  this  section  to  denote  the  reference  points  1  and  2 
since  the  region  index  q  is  not  needed.  This  assumes  that 
points  1  and  2  are  in  the  same  region.  Also,  quantities  with 
asteriks  in  this  section  represent  the  undebulked  configuration 
and  should  not  be  confused  with  the  previous  meaning  of  the 
asteriks. 
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and  (77).  A  good  first  approximation  is  z2  *  =  z2 .  The  constant 
c*  can  be  adjusted  for  a  closer  fit  between  the  initial  ply  boundary 
and  the  control  surface.  In  this  case,  a  revised  solution  of  (77), 
(79),  and  (80)  must  be  executed,  while  02*  is  given  by  (75). 

At  this  time,  we  have  the  coordinates  of  each  reference  point. 
We  also  have  the  mapped  image  of  each  reference  point  in  two 
coordinate  systems.  Thus,  the  next  step  is  to  establish  the 
relationship  between  the  two  coordinate  systems  such  that  the 
two  images  of  each  reference  point  will  coincide.  Referring  to 
Figure  7,  this  is  accomplished  by  use  of 

tan  0  sg*+TT*  (81) 

where 


S  =  R2  cos  ©2  -  Rj^  cos  0^^ 
T  =  R2  sin  ©2  -  Rj^  sin  ©^ 


(82) 


as  well  as 

x  =  R^  cos  0^  -  R^*  cos  (0^*  -  0)  =  R2  cos  ©2-R2*cos  ( O2*-0)  (83) 

and 

y  =  R1sin01-R1*sin(01*-0)=  R2sin02-R2*sin (  02*-  0)  (84) 

Given  an  arbitrary  image  (R,  0)  in  the  final  (cured)  state,  there¬ 
fore,  we  can  define  its  counterpart  (R*,  0*)  in  the  undebulked 
state  via 

R*  =  [(Rcos0-  x)  2  +  (Rsin0-  y)2]*5  (85) 


and 


tan  ( © 


i*  - 


0)  - 


Rsin©- 


3r 


(86) 


Rcos  ©- 

Finally,  we  arrive  at  the  undebulked  coordinates  (r*,  0*,  z*)  from 
the  relations 
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r.  , 

(A* )  +1 


(— )  1' 
'A*1  ‘ 


(87) 


z*  =  z^  + 


c*  [  (A*)2+11^ 
(A*)2 


(0*-^*  +  A*  +  cot  A*  -  Al*  -  cot  A  L* )  (88) 


9*  =  0  *  +  — j-  (z*-z.)  +  a,*  +  cota.*  -  a*  -  cota*  (89) 
l  c*  11  1 

which  are  derived  from  (8),  (9),  and  (2).  By  use  of  the  present  model, 

one  can  approximate  quite  accurately  the  condition  that  a  ply 
boundary  remains  in  contact  with  a  fabrication  tooling  surface, 
provided  the  ply  boundary  contains  no  discontinuity  in  the  region. 

The  above  treatment  defines  the  debulking  kinematics  of  an 
exact  involute  body.  For  a  body  constructed  by  the  start-line 
method,  eqs.  (75)  -  (89)  are  only  valid  within  a  strip  of  the  hypo¬ 
thetical  involute  surface.  Because  of  the  inherent  discontinuities 
in  the  hypothetical  involute  surface,  however,  the  form  of  the  most 
accurate  approximate  scheme  to  treat  the  debulking  kinematics  for  a 
body  generated  by  the  start-line  method  is  not  evident.  Therefore, 
these  problems  should  probably  be  examined  via  a  case  by  case 
analysis. 
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CONCLUDING  REMARKS 


Construction  of  bodies  such  as  rocket  nozzle  exit  cones  by 
the  involute  method  involves  analytical  modeling  of  various  aspects 
of  the  fabrication  procedure.  Furthermore,  the  fabrication  approach 
may  involve  the  use  of  an  exact  involute  ply  pattern  or  one  generated 
by  the  (inexact)  start-line  method,  which  may  include  interleaving. 
Since  the  start-line  model  contains  inherent  discontinuities,  it 
is  important  to  evaluate  the  effect  of  the  approximation  on  ply 
pattern  distortion  and  deformation.  Finally,  proper  debulking 
requires  carefully  designed  tooling  and  a  model  describing  the 
kinematics  of  the  debulk  process.  Aside  from  the  manufacturing 
considerations,  the  stress  analyst  needs  a  model  to  characterize 
the  variation  of  the  important  involute  parameters,  such  as  arc 
and  helical  angles,  throughout  the  body,  as  these  parameters 
define  the  distribution  of  the  elastic  stiffness  tensor,  which  in 
turn  influence  the  stress  field.  All  these  topics  have  been 
treated  here  including  the  development  of  explicit  relations  where 
practical.  Although  some  of  these  equations  appear  elsewhere, 
such  as  [1]  -  [4],  they  have  also  been  included  here  for  completeness 
and  convenience  of  the  reader. 
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Figure  7.  Transformation  for  Debulk  Analysis. 


APPENDIX  A 


A  NOTE  ON  SIGN  CONVENTIONS 


Since  many  parameters  are  involved  in  the  definition  of  the 
exact  involute  surface,  it  is  not  surprising  that  various  sign  con¬ 
ventions  appear  in  the  literature.  In  other  cases,  inconsistencies 
in  sign  are  present,  e.g.,  eq.  (B— 3)  of  [1]  and  eq,  (3)  of  [2]. 
These  have  been  corrected  in  the  corresponding  equations  of  the 
present  work.  Since  the  use  of  these  relations  to  define  stiff¬ 
ness,  stress,  and  displacement  components  in  stress  field  models 
requires  rigid  adherence  to  a  consistent  sign  convention,  we  shall 
review  the  present  nomenclature  in  explicit  fashion  at  this  point. 

All  equations  derived  in  this  work  are  consistent  with  the 
following  convention: 


a)  Unit  vectors  £,],£  lie  in  the  directions  of  (increasing)  0,r,z, 
respectively.  The  unit  vectors  i,j,k,  in  that  order,  form  a  right- 
hand  triad.  Furthermore,  the  involute  surface  is  assumed  to  be 

J  „ 

oriented  such  that  the  slope  of  the  meridian,  g^-,  is  not  negative. 

b)  Unit  vectors  t,s,n  must  be  oriented  such  that  they  are  defined 

by 


a 

t 

A 

S 


j  sin  y  + 


cos  y 


A  /V 

i  cos  a  -  j  sin  a 


n  sin  tp  =  t  x  s 

and  the  ranges  of  the  various  angles  are  given  by 

0  i  a  S  ir/2  ,  0  £  Y  £  ir/2  ,  tt/2  £  ip  i  it 


(A-l) 


(A-2) 
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A 

In  this  way,  n  is  the  unit  outward  normal  vector  of  the  involute 
surface . 

c)  The  helical  angle  $  may  be  chosen  to  represent  the  orientation 
of  any  particular  line  segment  on  the  ply  pattern.  For  definite¬ 
ness,  let  us  assume  $  is  measured  to  the  warp  direction.  Then  a 
positive  angle  <{>  is  measured  from  the  meridian  toward  the  warp 
direction  such  that  the  vector  representing  $  by  the  right-hand 
rule  lies  in  the  direction  of  n. 
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